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ABSTRACT 

We present the first evidence of clear signatures of tidal distortions in the density 
distribution of the fascinating open cluster NGC 6791. We used deep and wide-held 
data obtained with the Canada-France-Hawaii-Telescope covering a 2° x 2° area around 
the cluster. The two-dimensional density map obtained with the optimal matched filter 
technique shows a clear elongation and an irregular distribution starting from ~ 300" 
from the cluster center. At larger distances, two tails extending in opposite directions 
beyond the tidal radius are also visible. These features are aligned to both the absolute 
proper motion and to the Galactic center directions. Moreover, other overdensities 
appear to be stretched in a direction perpendicular to the Galactic plane. 

Accordingly to the behaviour observed in the density map, we find that both the 
surface brightness and the star count density profiles reveal a departure from a King 
model starting from ~ 600" from the center. 

These observational evidence suggest that NGC 6791 is currently experiencing 
mass loss likely due to gravitational shocking and interactions with the tidal field. We 
use this evidence to argue that NGC 6791 should have lost a significant fraction of 
its original mass. A larger initial mass would in fact explain why the cluster survived 
so long. Using available recipes based on analytic studies and IV-body simulations, 
we derived the expected mass loss due to stellar evolution and tidal interactions and 
estimated the initial cluster mass to be M- ln i = (1.5-4) x 1O 5 M 0 . 
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1 INTRODUCTION 

NGC 6791 is one of the most intriguing open cluster (OC) 
in the Galaxy. It possesses a unique set of properties that 
have attracted a lot of attention in the last decade. 

Indeed, NGC 6791 is one of the most massive (M ~ 
5000 JMq; Platais et al 2011), oldest (t ~ 8 Gyr; King et al. 
2005; Grundahl et al. 2008; Garcia-Berro et al. 2010; Bro- 
gaard et al. 2012) and most metal-rich ([Fe/H] ~ +0.40; 
Carraro et al. 2006, Gratton et al. 2006, Origlia et al. 2006) 
Galactic OCs. More specifically, this combination of prop¬ 
erties is unique for clusters at its Galactocentric distance 
(d ~ 8A'pc; Carraro 2014). Stars with similar age and metal- 
licity are found in the Galactic bulge (Bensby et al. 2013) 
where it has been suggested (Jilkova et al. 2012) NGC 6791 
might have formed. 

Given its peculiarities, and in particular its mass, 
NGC 6791 represents an ideal link between OCs and glob¬ 


ular clusters. For this reason it has been widely targeted to 
understand whether star-to-star variations of light element 
(C, N, O, Na, O, Mg and Al) abundances can be found also 
in this system as observed in basically all globular clusters 
(Gratton et al. 2012). This topic is still matter of debate and 
no consensus has been reached so far on whether the clus¬ 
ter hosts more than one generation of stars. Geisler et al. 
(2012) presented evidence of Na-O anti-correlation among 
giant stars in NGC 6791, while recently Bragaglia et al. 
(2014) and Cunha et al. (2015) have questioned this result 
finding an homogeneous Na abundance (within ~ 0.1 dex). 

NGC 6791 shows also an anomalous horizontal branch 
(HB) with a well populated red clump, as expected for its 
high metallicity, and a group of hot HB stars (Liebert et al. 
1994, Tofflemire et al. 2014), whose origin is still not com¬ 
pletely understood (Carraro 2014). The combined UV flux 
of the few hot HB stars and its metallicity makes this cluster 
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a fairly good proxy of standard UV-upturn elliptical galax¬ 
ies (Buzzoni et al. 2012). This result is supported also by 
the Lick indices analysis of its spectral energy distribution. 

This plethora of unique properties makes NGC 6791 
an extremely interesting object to study and understand. 
A number of questions about the nature of this system 
still arise: how and where could such a stellar system have 
formed? Is NGC 6791 a genuine OC? Did it form close to the 
bulge? How could have survived for so long in the adverse 
high-density environment of the inner Galactic disk, where 
it presently is? One of the still missing key ingredient that 
could help to understand the origin of such a fascinating 
system is its original mass. 

The mass budget of any stellar system is strongly af¬ 
fected by two-body relaxation, stellar evolution, encounters 
with giant molecular clouds and spiral arms, interactions 
with the tidal field. The disruption scenarios and the sur¬ 
vival rates of Galactic OCs have been studied since the 
sixties (Spitzer & Harm 1958; Spitzer 1958). Recently, N- 
body simulations (e.g., Gieles et al. 2006; Gieles & Baum- 
gardt 2008) have also shown the impact of giant molecular 
clouds and of the initial structure of OCs on their dissolu¬ 
tion. The OC orbits have generally small eccentricities and 
they tend to be located at a small height ( Z ) on the Galac¬ 
tic plane, as a consequence they are expected to pass many 
times trough the Galactic disk. Each of these gravitational 
shocks heat up and compress the cluster which then takes 
an elongate shape which flattens to its maximum at Z = 0 
(Leon 1998). Repeated disk-shocking speed-up the disrup¬ 
tion of OCs (Combes et al. 1999) and after each passage 
cluster stars rejected in the halo of the system are stripped 
out by the gravitational field of the Galaxy. These ejected 
stars are expected to form tidal tails extending far from the 
cluster inner regions. Indeed, Davenport & Sandquist (2010) 
and Bergond et al. (2001) have found that some OCs show 
a quite elongated and extended structures likely resulting 
from the interactions with the Galactic tidal field. 

In general, OCs are thought to dissipate in a timescale 
of the order of ~ 10 7 — 10 8 yr (Binney & Tremaine 1987). 
In this context the existence of very old OCs like NGC 6791 
represents an anomaljQ. This system may have managed to 
survive so long either because it moved along a preferential 
orbit or its original mass was much larger than observed 
today. 

In this paper we analyse the two-dimensional (2D) and 
projected density distributions of NGC 6791 with the aim 
of understanding whether the cluster suffered from strong 
tidal interactions that could have significantly decreased its 
mass. We find that this system has an elongated structure 
and shows overdensities in the external regions that are 
aligned to the direction of the orbital motion. These fea¬ 
tures strongly suggest that this system has likely lost part 
of its mass because of gravitational shocks and interaction 
with the Galactic tidal held. We use this argument to sup¬ 
port the hypothesis that this system was significantly more 
massive than observed today and we derive an estimate for 
its original mass by means of a simplified analytic approach. 

The paper is organised as follows: in Section 2 data-set 


1 Other old (t > 2 — 3 Gyr) OCs for which such an argument 
applies are NGC 188, M 67 and Berkley 11. 
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Figure 1. Map of the entire CFHT database with respect to the 
position of C gra v (Xq.Yq). The two circles represent the location 
of core (r c ) and tidal radii (rt). The box in the upper right corner 
defines the "control field” (CF) area. 

and data-reduction procedures are described, in Section 3 
the 2D density map is obtained while in Section 4 both the 
star-count projected density and the surface brightness pro¬ 
files are analysed; finally in Section 5 we estimate, the initial 
mass of NGC 6791. 


2 OBSERVATIONS AND DATA REDUCTION. 

We used deep proprietary images (Prop ID: OPTICON 
2013B/005; PI Moitinho) obtained with the wide field im¬ 
ager MegaCam mounted at the Canada-France-Hawaii Tele¬ 
scope (CFHT). This camera consists of a mosaic of 36 chips 
of 2048x4612 pixels each, with a pixel scale of 0.185” pixel -1 
providing a total field of view of ~ 1° x 1°. We used four 
different pointings set up to cover a total field of view of 
~ 2° x 2° centered approximately on the cluster center (Fig¬ 
ure 1). For each pointing four images acquired both with the 
g' and r' bands and with exposure times t exp = 100 sec and 
texp = 180 sec respectively have been secured. An appro¬ 
priate dither pattern of few arcesc have been adopted for 
each pointing. This allowed us to cover most of the inter¬ 
chip gaps, with the exception of the most prominent and 
horizontal ones (see Figure 1). 

The data were pre-processed (i.e. bias and flat-field cor¬ 
rected) by the Elixir pipeline developed by the CFHT team. 
By means of an iterative procedure, an adequate number 
(> 20) of isolated and bright stars have been selected in 
each chip and band to model the Point Spread Function 
(PSF). Then the PSF model was applied to all the stellar- 
like sources at about 4<r from the local background by us¬ 
ing DAOPHOT and the PSF-htting algorithm ALLSTAR 
(Stetson 1987). For each hlter and chip we matched the 
single-frame catalogues to obtain a master list. Finally, the 







Tidal distortions around the old open cluster NGC 6791 3 


g' and r' master lists were then geometrically matched and 
combined. Each master list includes the instrumental mag¬ 
nitude, defined as the weighted mean of the single image 
measurement reported to the system reference frame of the 
transformation, and the error, which is the standard error 
of the mean. 

The instrumental coordinates have been reported inde¬ 
pendently for each chip to the absolute astrometric system 
(a, S ) using a large number of stars in common with the 
DRY release of the Sloan Digital Sky Survey (Abazajian et 
al. 2007) and following the procedure described by Dalessan- 
dro et al. (2009). Using the same stars, we derived for each 
chip and band the zero-points to report the instrumental 
magnitudes to the Sloan photometric system. At this stage 
all the master lists are in the same photometric and astro- 
metric system. We then merged them to create the final 
catalogue, which counts about 270,000 stars. 

The ( gg' — r') colour-magnitude diagrams (CMDs) of 
NGC 6791 and of the surrounding field are shown in Fig¬ 
ure 2. In the CMDs, it is possible to clearly distinguish the 
Main Sequence (MS) of the cluster (left panel of Figure 2), 
extending down to g' ~ 23.5 in the innermost regions, cor¬ 
responding to about five magnitudes below the Turn-Off. 
Also, it is possible to note that our data suffer of strong sat¬ 
uration starting at g' ~ 17. Two almost vertical sequences 
located around (</ — r') ~ 0.6 and (</ — r') ~ 1.5 are clearly 
visible both in the innermost (Figure 2 left panel) and ex¬ 
ternal regions (Figure 2 right panel). They are populated by 
M-dwarf and Main-Sequence Galactic disk stars. 

A differential reddening of about A E(B — V) ~ 0.05 has 
been estimated for this cluster (see for example Platais et al. 
2011). Indeed, at the level of the Turn-Off (g 1 ~ 18), the MS 
shows a mild broadening not explainable in terms of pure 
photometric errors. Such an effect disappears at fainter mag¬ 
nitudes (</ ~ 19), where the MS starts to bend diagonally 
in the (g 1 , g' — r') CMD. Given its relatively small ampli¬ 
tude, we will neglect the effect of differential reddening in 
the following analysis. 

3 THE DENSITY MAP 

3.1 The matched filter analysis 

To probe the spatial distribution of NGC 6791 we first anal¬ 
ysed the 2D density distribution. 

The density analysis of NGC 6791 has been performed 
in the circular area enclosed within a distance from the clus¬ 
ter center (see details in Section 4) r = 3000”, thus covering 
a 1.7° x 1.7° area. We selected stars along the MS in the 
magnitude range 18 < g' < 21.5 and within three times 
the local color dispersion about the cluster mean ridge line 
(see Figure 2). Stars selected in this way define the ’’Cluster 
Population” (CP) sample. These selection criteria have been 
chosen to ensure a good level of completeness and to mini¬ 
mize the contamination from Galactic field stars. Variations 
of the adopted magnitude limits do no affect the final result. 

As a first qualitative step, we obtained a rough 
smoothed density distribution of CP stars. As expected, at 
a first look, the distribution of CP stars appears to be ex¬ 
tremely concentrated around the position of the cluster and 
it shows some irregularities and structures allover the field 
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Figure 2. (</, g'—r 1 ) CMDs of the innermost region of NGC 6791 
(left panel) and of the ’’Control Field” (right panel). The black 
box along the MS of NGC 6791 encloses the region in the CMD 
where CP stars have been selected. The shaded area defines the 
region where fiducial Galactic disk MS stars have been selected. 

of view (FOV). We also note that the North-West quadrant 
shows on average a smoother distribution than the other 
quadrants for r > 1000”. 

To investigate in more detail the two-dimensional dis¬ 
tribution of CP stars, we adopted the widely used opti¬ 
mal matched filter technique (see for example Kuhn et al. 
1996; Rockosi et al. 2002; Odenkirchen et al. 2001, 2003). 
It has been shown (Davenport & Sandquist 2010) that this 
method is robust for identifying low-density stellar popu¬ 
lations against a significant background, as in the case of 
OCs. The optimal matched filter technique requires to use 
a fiducial field population as reference. For this purpose, 
we selected stars in the same magnitude and colour bins as 
the CP sample and located at r > 3000” from the cluster 
center in the North-West quadrant (because of its apparent 
homogeneous distribution). We named this sample of stars 
“Control Field” (CF; Figures 1 and 2 ). 

We computed the densities in the CMDs (Hess dia¬ 
grams) of both the CP and CF samples. Then we assigned 
to each star in the CP sample a weight defined as the recip¬ 
rocal of the fraction of CF stars lying within a fixed range 
of colour and magnitude from its location in the CMD. 

Finally, the weighted distribution of star positions 
was transformed into a smoothed surface density function 
through the use of a kernel whose width has been chosen 
to be large enough so as to make negligible the impact of 
the gaps in the FOV. This procedure yields the surface den¬ 
sity distribution shown in Figure 3. As apparent, the overall 
density distribution is elongated and irregular. In particular, 
we observe that the innermost ~ 300” are quite spherical, 
while at larger distance the cluster starts to be distorted 
mainly along the North-East/South-West direction in an al¬ 
most symmetrical way respect to the center of the cluster. 
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Figure 3. Large-scale two-dimensional colour-coded surface density map around NGC 6791. The contour levels span from 3c r to 40cr 
with irregular steps. The solid arrow represents the direction of the absolute proper motion of NGC 6791, while the dashed ones mark 
the direction of the Galactic center (GalC) and that perpendicular to the Galactic plane (Z=0). 


At r > 1000”, two tails running at opposite directions are 
clearly visible, with the one in the southern direction extend¬ 
ing up to distances larger than 2500” from the cluster center. 
It is important to note that both the external tidal arms and 
the innermost elongation of the cluster run almost parallel 
to the absolute proper motion vector (Bedin et al. 2006) and 
to the direction of the Galactic center (GalC; Figure 3). This 
behaviour is in agreement with a real tidal nature for the ob¬ 
served distortions. In fact, the orbital motion of NGC 6791 
likely develops close to the Galactic equatorial plane (Jflkova 
et al. 2012) and, since tidal tails are supposed to lie on the 
orbital plane (at least in axisymmetric potentials, see Mon- 
tuori et al. 2007), they are expected to be seen with such a 
projected alignment. Moreover, other distortions extending 
almost perpendicularly to the PM direction are clearly vis¬ 
ible and their formation could be related to the tidal effect 
due to the Galactic disk (e.g., Bergond et al. 2001). 

The features observed in the 2D density distribution of 
NGC 6791 are indeed typical evidence of recent stellar mass 
loss due to dynamical evolution of the system under the tidal 
influence of the overall Galactic potential. It is interesting 
to note that the density map of NGC 6791 is qualitatively 
similar to that observed in M 67 by Davenport & Sandquist 
(2010) and to a weaker extent to other three OCs (namely 
NGC 2287, NGC 2516 and NGC 2548) by Bergond et al. 
( 2001 ). 

We also note that in the South-East quadrant for 
r > 2000” there is a star density increase that does not 


seem to be directly connected to the cluster. Interestingly 
enough, this overdensity appears at decreasing Galactic lat¬ 
itude therefore we speculate that it might be due to the 
density gradient of Galactic disk stars. 

3.2 Sanity checks 

To further check the significance of the anisotropic dis¬ 
tribution of NGC 6791 we used for comparison stars se¬ 
lected in the CMD in the range 0.4 < ( g' — r') < 0.8 and 
20 < g' < 21.5 (see Figure 2). These are likely MS Galactic 
disk stars and they trace the behaviour of the main contrib¬ 
utors to the contamination of CP stars. This sample is sup¬ 
posed to be homogeneously distributed across the FOV and 
actual deviations from a smooth distribution reflect varia¬ 
tions in the detection efficiency, Galactic field and/or extinc¬ 
tion gradients. 

We performed on this sub-sample the same analysis 
done previously for CP stars. The 2D density map of MS 
Galactic disk stars is shown in Figure 4. We observe that 
in general they are uniformly distributed across the FOV 
as expected, but they show a density excess in the external 
regions of the South-East quadrant. This evidence further 
suggests that the relevant feature in the 2D density distri¬ 
bution of CP stars at r > 2000” (Figure 3) in this quadrant, 
is likely due to a gradient of the density of the Galactic disk 
within the FOV. 

In addition, we divided the FOV in ten angular sec- 
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Figure 4. As in Figure 3, but for stars with 0.4 < (g' — r') < 0.8 
and 20 < g' < 21.5, which are representative of the Galactic field 
contamination (see Section 3.2). Colors are in the same relative 
scale as in Figure 3. 

tors of 36° each with vertex at the cluster center. Then for 
each sector with angle 6 measured counter-clockwise from 
the West direction, and for different distances from the clus¬ 
ter center we estimated the ratio between CP stars ((Vcp), 
and likely Galactic disk MS stars, selected as described be¬ 
fore (IVfieid). The results for r < 500" and 500" < r < 1500" 
are shown in Figure 5. In both cases the ratio shows sig¬ 
nificant variations as a function of 8, as expected for a real 
anisotropic density distribution. Moreover it is also worth 
noting that different patterns in the central and outermost 
regions appear. For 500" < r < 1500" the maxima of the 
ratio occur at 9 = 100° and 9 — 280° corresponding to the 
two opposite main tidal arms. 

These sanity checks demonstrate that the observed 2D 
density distribution shown in Figure 3 is a real feature of 
NGC 6791 and cannot be ascribed to observational bias or 
field gradients. 


4 STAR DENSITY AND SURFACE 
BRIGHTNESS PROFILES. 

Star density distortions, elongation and tidal tails are ex¬ 
pected to be detected also in the radial density profile (see 
for example Odenkirchen et al. 2003; Capuzzo-Dolcetta et 
al. 2005). Therefore as a second step of our analysis, we 
studied the projected density profile of NGC 6791. 

First we derived the center of gravity (C gra v) by using 
an iterative procedure (see for example Dalessandro et al. 
2013) and averaging the positions a and 5 of CP stars. We 
used the center reported by the WEBDA web pag43 as start- 

2 http://www.univie.ac.at/webda/cgi- 
bin/ocl_page.cgi?dirname=ngc6791 


Figure 5. Ratio between CP and field stars as a function of the 
position angle 9 (Section 3) for stars at r < 500" (upper panel) 
and at 500" < r < 1500" (lower panel). 

ing guess point of our procedure. To avoid spurious effects 
due to incompleteness we performed several estimates by 
adopting different magnitude and distance selections. In par¬ 
ticular, we used three different magnitude bins in the range 
18 < g' < 21.5, with the upper limit (fainter magnitude) 
increasing by 0.5 mag at each step. For each magnitude bin, 
we repeated the procedure for four different radial selections, 
from 250" to 400" with a step of 50". We ended up with a 
total of 20 estimates. The resulting C grav is the average of 
these measures and it is located at a c = 19 h 20 m 54 s .273 and 
S c = 37°46'25".20 (R.A.= 290.2261395 Dec=37.77366712), 
with an uncertainty of ~ 1.2". This new determination is 
located at ~ 12" North-East from the center reported by 
the WEBDA web page. While it is not trivial to understand 
the origin of this discrepancy, we stress here that it has a 
negligible impact on the results presented in this paper. 

The projected density profile of NGC 6791 has been de¬ 
termined using direct counts of CP stars. We excluded the 
South-East quadrant from the analysis since in the most ex¬ 
ternal part of this region the star distribution is dominated 
by MS Galactic disk stars (Section 3.1). Using the proce¬ 
dure described in Dalessandro et al. (2013), we divided the 
considered FOV into 18 concentric annuli, each centered on 
C gr av and suitably split in an adequate number of subsec¬ 
tors. Number counts have been calculated in each subsec¬ 
tor and the corresponding densities were obtained dividing 
them by the sampled area. Incomplete area coverage, mainly 
due to the presence of gaps (Figure 1), have been properly 
taken into account. The stellar density of each annulus was 
then defined as the average of the subsector densities and its 
standard deviation computed from the variance among the 
subsectors. The resulting projected surface density profile is 
shown in Figure 6. 

We estimated the background density contribution by 
using only stars in the North-West quadrant. In fact, as 
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discussed before (see also Figure 3) this region is free by 
any evident structure caused either by the cluster elonga¬ 
tion or from Held stars. The background contribution has 
been estimated as the average of the four outermost density 
measurements (r > 1400”), which define a kind of plateau 
at log(pbck) ~ —3.6 stars arcsec -2 . We then subtracted this 
value to the observed density profile. We have verified that 
this choice is appropriate. In fact, the contribution to the 
background due to field stars (see selection criteria in Sec¬ 
tion 3.1) is constant in all quadrants but the South-East 
one. 

The derived density profile has been reproduced by us¬ 
ing a single-mass King model (King 1966). The fitting pro¬ 
cedure is fully described in Miocchi et al. (2013). The fit 
was limited to the innermost r < 500” to avoid the most 
elongated regions of the cluster. The best-fit is obtained for 
a concentration c = 1.06 (corresponding to a central di¬ 
mensionless potential Wo = 5.3) and core radius r c = 160”. 
These parameters are in partial disagreement with those ob¬ 
tained by Platais et al. (2011; c = 0.74, r c = 196”) which 
to our knowledge is the only other density profile analy¬ 
sis present in the literature for NGC 6791. The authors 
obtained the density profile using stars with g' < 22 and 
weighted according to their proper motion. However their 
data were limited to the innermost 900” and therefore the 
determination of the tidal radius and any residual in the 
background contribution was rather unstable. 

As apparent in Figure 6, for log(r) > 2.8 (r > 600”) the 
density profile clearly deviates from the behaviour predicted 
by the King model. The use of a Wilson model (Wilson 1975) 
does not provide any significant improvement on the overall 
quality of the fit. We have also checked that even if we do 
not constrain the model to fit the innermost region, there 
is no way to satisfactory reproduce the density profile for 
r > 600”. The external part of the projected density pro¬ 
file instead follows a power-law behaviour with an exponent 
a ~ —1-7, in agreement with the values found in Galactic 
globular clusters with observed tidal tails (see for example 
Sollima et al. 2011 and references therein) and as predicted 
by theoretical models (Johnston et al. 1999). Kfipper et al. 
(2010) analysed the variation of the power-law slopes of tidal 
tails of dissolving clusters as a function of the orbital phase, 
eccentricity and initial cluster densities by means of N-body 
simulations. They found that the power-law slopes have typ¬ 
ical values a = —4-5 that can increase reaching values 

a ~ —2 for systems with highly eccentric orbits and close to 
or at the apogalacticon, where tidal tails get strongly com¬ 
pressed. This upper limit is in good agreement with what 
found for NGC 6791, however we should note that the anal¬ 
ysis by Kfipper et al. (2010) is performed on clusters moving 
on orbits perpendicular to the line of sight, while this is not 
the case for NGC 6791. As a consequence the relatively shal¬ 
low power-law observed in NGC 6791 can be the result also 
of projection effects. 

To further investigate the results just described and 
check whether they might be caused by incompleteness or 
bias of our photometric catalogue, we also derived the sur¬ 
face brightness profile by performing the analysis directly 
on FITS images (see for example Dalessandro et al. 2012). 
We used a low-resolution single-plate ESO Digitalized Sky 
Survey image obtained in optical band and extending up 
to about 2000” from C gra v We limited the analysis to the 



Figure 6. Observed star count density profile as a function of ra¬ 
dius (open grey squares). The dashed line represents the density 
value of the background as obtained by averaging the four out¬ 
ermost density measurements in the North-West quadrant. The 
black filled dots are densities obtained after background subtrac¬ 
tion. The best single-mass King model is also overplotted to the 
observations (solid line) and the structural parameters are la¬ 
beled. The lower panel shows the residual between the observa¬ 
tions and the best-fit model. 


same quadrants used for deriving the star counts density 
profile and we used the North-West quadrant to estimate 
the background contribution. We divided the FOV in a sim¬ 
ilar number of concentric annuli and subsectors. The surface 
brightness assigned to each annulus is given by the ratio be¬ 
tween the mean of the fluxes measured in the subsectors 
and the area covered by the annulus. The resulting surface 
brightness profile is shown in Figure 7. First we observe that 
the King model best-fitting the star counts density profile 
well reproduces also the surface brightness profile. More im¬ 
portantly, we note that consistently with what found in the 
star count projected density profile, for r > 600” the sur¬ 
face brightness profile starts to deviate from the behaviour 
predicted by the model and it declines as a power-law. 

We can therefore safely conclude that both the 2D den¬ 
sity map and the density (surface brightness) profile analysis 
suggest that NGC 6791 is currently experiencing tidal strip¬ 
ping events. 

It is worth noting that King et al. (2005) observed a 
rather flat mass function in the centre of NGC 6791 by using 
very deep Hubble Space Telescope observations. This obser¬ 
vational evidence further supports the results obtained in 
this work. In fact, for clusters that experienced episodes of 
mass loss, a flat mass function is expected because of the 
preferential depletion of low mass stars. A detailed analysis 
of the radial variation of the mass function of this system 
will be presented in an upcoming paper (E. Dalessandro et 
al. 2015). 
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Figure 7. Optical instrumental surface brightness profile for 
NGC 6791.The grey empty squares represent the observed raw 
profile, while the black circles are the sky subtracted values. The 
horizontal dashed line marks the estimated of the background 
level as inferred in the North-West quadrant. The background- 
subtracted profile is shown with filled black circles. 


5 MASS LOSS AND INITIAL MASS 

ESTIMATE. 

The observational evidence presented in the previous Sec¬ 
tions reveal that NGC 6791 might have lost a fraction of its 
original mass due to environmental effects (like tidal shocks 
due to close interactions with giant molecular clouds, spiral 
arms, the Galactic disk and, in general, to the interactions 
with the Galactic tidal field) or internal dynamical effects 
(like two-body relaxation. In this Section we attempt to es¬ 
timate the total mass lost by NGC 6791 during its evolution 
and derive its original mass. 

While many recipes can be used for this purpose (e.g., 
Vesperini et al. 2013), we adopted the approach described 
by Lamers et al. (2005). It has the advantage of describing 
the way the mass of a cluster decreases with time by means 
of relatively simple analytic expressions. The analysis by 
Lamers et al. (2005) accounts for the effect of both stellar 
and dynamical evolution. It is important to note that the 
results obtained with this simplified approach are in good 
agreement with those obtained by detailed TV-body simula¬ 
tions for clusters in the tidal field of the Galaxy (see also 
Lamers & Gieles 2006). 

Lamers et al. (2005) provide an expression for the initial 
mass of the cluster, Mini, in terms of its age (t) and present 
mass (M): 


Mini ~ 


M 

M 0 


7 i 

to 


n 1/7 


[1 Qev(t)] 


(1) 


to is the dissolution time scale parameter. It is a constant 


depending on the strength of the tidal field and it basically 
describes the mass loss due to dissolution processes such as 
Galactic tidal field interactions and shocks due to encoun¬ 
ters with giant molecular clouds or spiral arms (Gieles et 
al. 2006, 2007). The smaller values of to are typically asso¬ 
ciated to encounters with molecular clouds and spiral arms, 
while the larger to the Galactic tidal field. Also, the longer 
to, the weaker the tidal effect is. We used to = 3.31 ^'q Myr 
(corresponding to the disruption time of 1.3 ± 0.5 Gyr for a 
cluster with the initial mass of 10 4 Mq), obtained by Lamers 
et al. (2005) by comparing the distribution of mass and age 
of OCs in the solar neighbourhood with theoretical predic¬ 
tions. This choice is reasonable since the current Galacto- 
centric distance of NGC 6791 roughly corresponds to that 
of the Sun. 

7 is a dimensionless index, which depends on the clus¬ 
ter initial density distribution and has typical values ranging 
from about 0.6 to 0 . 8 , which increases with the concentra¬ 
tion, as constrained from theoretical studies (Gieles et al. 
2004) and observations (Boutloukos et al. 2003). 7 is usually 
constrained by the King dimensionless potential Wo, which 
characterises the concentration of the King models (King 
1966). We adopted here 7 = 0.62 corresponding to Wo = 5 
that is a typical value for open clusters and is in agreement 
with what derived by the analysis of the observed star count 
density and surface brightness profiles (Wo = 5.3 ± 1.0; Sec¬ 
tion 4). However, it is important to note that the star cluster 
density distributions change with time, therefore including 
the present-day value of Wo in Eq. © represents a crude 
approximation. I 11 particular, it is well known (e.g. Trenti 
et al. 2010 ) that the concentration of stellar systems tend 
to increase (and Wo accordingly increases) as the dynamical 
evolution of star cluster proceeds. As a consequence, given 
the proportionality linking 7 to Wo, if the adopted value of 
Wo is larger that its initial one, also the adopted value of 7 
will represent an upper limit to the correct one. 

The function q ev ( t ) describes the mass loss due to stellar 
evolution and can be approximated by the following analyt¬ 
ical formula: 

log 10 q ev (t) = (log 10 t - a) b + c, for t > 12.5 Myr, (2) 

where a, b, and c are coefficients that depend on the cluster 
metallicity. The [Fe/H] abundance of NGC 6791 is about 
+0.4 (Z ~ 0.05; Carraro et al. 2006, Gratton et al. 2006, 
Origlia et al. 2006) leading (see Table 1 and Eq. 2 by Lamers 
et al. 2005) to values of a = 7.00, b = 0.25, and c = —1.82. 

Eq. Q depends also on the cluster current mass M and 
its current age t. We adopted a current mass M = 5000 Mg 
as derived by Platais et al. (2011) and age t = 8 Gyr (Grun- 
dahl et al. 2008, Garcia-Berro et al. 2010, Brogaard et al. 
2012). These values posses significant errors and their exact 
estimate is still a matter of debate, however it is impor¬ 
tant to note that the final error in the initial mass esti¬ 
mate is driven mostly by the uncertainty in the dissolution 
time scale parameter to- The dependence of the cluster ini¬ 
tial mass on to within the range of 3.31 i'q Myr is plotted in 
Fig- [SI together with the spread caused by the uncertainty of 
the current cluster mass (lower and upper limit of 4000 and 
6000 Mg considered) and its current age (lower and upper 
limit of 7 and 9 Gyr considered). 

The initial mass of NGC 6791 obtained by Eq. 0 turns 
out to be Mini = (1.5-4) x 10 5 Mq, i.e. more than 50 times 
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Figure 8. Cluster initial mass given by Eq. {TJ as a function of 
the dissolution time scale parameter to- The red curve shows the 
dependence derived for the current cluster mass of 5000 Mq and 
its age of 8 Gyr; the dark-grey coloured area corresponds to values 
derived for mass within 4000-6000 Mq at age of 8 Gyr; the light- 
grey area corresponds to values derived for the mass of 5000 Mq 
at an age within 7-9 Gyr. 

larger than its present-day mass. For reasonable smaller val¬ 
ues of wE the initial mass derived by means of Eq. © will 
be slightly larger, but still compatible within the errors with 
the value quoted above. 

It is important to stress that this estimate is based on 
a simplified approach and on parameters derived by the av¬ 
erage behaviours of OCs. Different recipes and assumptions 
may likely lead to slightly different results. 


6 CONCLUSIONS 

In this study we have presented a deep, wide field, photo¬ 
metric data-set extending well beyond the tidal radius of the 
massive open cluster NGC 6791. This data-set allowed us to 
derive updated estimates of its structural parameters, by 
fitting a King profile over both the star counts and surface 
brightness profiles. With respect to previous studies, that 
were limited in spatial coverage (Platais et al. 2011), the 
King model best fitting the projected density profile turns 
out to have larger values of both c and rt- 

More interestingly, for the first time, we detected clear 
signatures of tidal features in the star distribution of 
NGC 6791, in the form of irregular, but evident, elonga¬ 
tion and tidal tails. Some of them seem to closely follow 
the cluster motion, while others are nearly perpendicular 
to the Galactic disk. These are clearly visible both in the 
2D surface density map and in the projected density profile 
and they represent typical evidence of recent stellar mass 
loss. We can therefore argue that NGC 6791 might have lost 

3 Note that OCs with very small values of the dimensionless po¬ 
tential (Wo < 3), are expected to be completely destroyed in a 
short timescale (~ 10 7 yr). 


quite a significant fraction of its original mass while orbit¬ 
ing around the Galactic center. By using a simple analytic 
approach (Lamers et al. 2005), we estimated the mass likely 
lost by NGC 6791 during its evolution because of the ef¬ 
fect of both stellar evolution and dynamical interactions. 
From this calculation we estimated the mass at its birth to 
be Mini = (1.5-4) x 10 5 Mq, several tens larger than its 
present-day mass. 

This derivation lends support to a scenario in which 
NGC 6791 was a much more massive system. This would 
qualitatively explain why the cluster could have survived 
for such a long time (its age is around 7-8 Gyr) contrary 
to the expectations of current estimates for the destruction 
rates of OCs in the Galaxy. 

It is interesting to note also that Mini is comparable 
to the present-day mass of most Galactic globular clusters 
and it is in rough agreement with the initial mass of the 
least massive clusters hosting multiple stellar populations 
(see for example Dalessandro et al. 2014). Therefore, the 
estimate of Mini obtained in this paper might set interesting 
constraints on the initial conditions of stellar systems able 
to form multiple stellar populations. 
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